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o 

We consider the projected SO(5) bosonic model introduced in order to connect the SO(5) the- 
ory of high-T c superconductivity with the physics of the Mott-insulating gap, and derive the cor- 
responding effective functional describing low-energy degrees of freedom. At the antiferromagnetic- 
^"1, superconducting transition, SO(5) symmetry-breaking effects due to the gap are purely quantum 

mechanical and become irrelevant in the neighborhood of a possible finite-temperature multicritical 
point separating the normal from the antiferromagnetic and the superconducting phases. A differ- 
ence in the magnon and hole-pair mobility always takes the system away from the 5*0(5) -symmetric 
fixed point towards a region of instability, and the phase transition between the normal and the two 
ordered phases becomes first order before merging into the antiferromagnetic-superconducting line. 
Quantum fluctuations at intermediate temperatures, while introducing symmetry-breaking terms in 
the case of equal mobilities, tend to cancel the symmetry-breaking effects in the case of different 
mobilities. 



o 

(N 



(N 



30 ■ IWC'S numbers : 7 !. 20. -z Ec 

c3 



X 



I. INTRODUCTION 

The 50(5) theory of high-T c superconductivityEHjl has been introduced as a concept to unify antiferromagnetism 
(AF) and d-wave superconductivity (SC) under a common symmetry principle. In order to study the physical 
i— 1 1 consequences, and to make predictions to compare with experiments, several exact 5*0(5) -symmetric models with 
small symmetry-breaking terms have been proposed and investigated in detailETQ. However, a shortcoming of these 
models, and of an exact 50(5) theory in general is that they are inconsistent with the antiferromagnetic gap at 
ly-^ , half filling, one of the most important features of the high-T c cuprates&tl. This can be understood by the fact that 
00 ' an 5*0(5) transformation "rotates" spin into charge and, thus, a requirement for an exactly 50(5) -invariant system 
would be to have the same charge and spin gap. This is in contradiction with the experimental situation in the high-T c 
I/"") ■ materials, where a large charge gap of some eV is present in the AE-statepat half-filling, while spin- wave excitations are 
ungapped. The introduction of a small symmetry-breaking termo'Enj'tjO, while on the one hand correctly selecting 
the AF state at half-filling and shifting the AF-SC transition to finite doping, does not introduce a charge gap of the 
correct order of magnitude. In contrast, in a weakly-coupled Hubbard ladder model, a spin and a charge gap q f t he 
same size are present and, in fact, it has been shown that 50(5) symmetry is dynamically restored at half fillingEHlEj. 

In order to cure this problem at strong coupling as well, a class of SO (5) models - "projected" 50(5) models - 
has been introduced where the Mott-Hubhard gap is taken into account by means of a Gutzwiller projection, whereby 
*^ ■ doubly-occupied states are projected outo. In that , it was shown that, despite the symmetry-breaking effects of the 
projection, static correlation functions remain exactly 50(5) symmetric within a mean-field approximation. This is 
O ■ due to the fact that, neglecting dynamic effects, the Hamiltonian is manifestly 50(5) invariant. However, dynamic 
effects breaki ng i th e 50(5) symmetry become important whenever quantum fluctuations are taken into accounto. In 
another paperl 15 H ll n, it was shown that the projection is crucial in order to correctly relate the <i-wave superconducting 
gap at finite doping with the <i-wave modulation of the AF gap observed at half filling by ARPES experimental^. 

In a microscopic physical system, one would in general expect 50(5) symmetry to be explicitly broken by several 
terms. This is certainly the case for the Hubbard model, for example. However, it often occurs in nature that a 
symmetry, which is broken on the microscopic level, is then restored in the long- wavelength limit. Concerning 50(5) 
, this has been shown to happen for quite generic ladder systems of the Hubbard typoljO. Recently, Murakami and 
Nagaosa argued that the bicritical point of the AF to SC transition in the organic superconductorpt^BEDT-TTF^X 
(Bis(ethylenedithio)tetrathiafulvakne) with X=Cu[N(CN)2]Cl, shows 50(5) critical exponental3. In fact, one of 
the scenarios suggested by ZhangEl is that there might be a direct first-order AF to SC transition terminating at a 
finite-temperature bicritical point, where the 50(5) symmetry is asymptotically restored at long wavelengthso. These 
ideas are very interesting from an experimental point of. view, and open the possibility of an explicit test of 50(5) 
symmetry, via a direct "measurement of the number 5"t2l. This could be done, as suggested in Ref. [l^, by measuring 
the critical exponents of the AF to SC transition, which, given the spatial dimensionality, should only depend on the 
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number of components n of the order parameter. On the other hand, it is toII known that for n > n c w 4, the 5*0(5) 
symmetric fixed point is unstable towards a so-called biconical fixed pointEj. However, since n = 5 is close to n c , it 
turns out that the stable biconical fixed point only breaks the symmetry by about 20%. 

However, the situation of the high-T c materials is quite delicate. As discussed above, the Mott-Hubbard gap plays 
an important role, and it produces a substantial breaking of the 50(5) symmetry. In Ref. it was shown that in the 
extreme case of a Gutzwiller projection, a degree of freedom is eliminated completely and the real and imaginary part 
of the local superconducting parameter become conjugate variables. Therefore, it is not clear whether such a projected 
50(5) symmetry can become asymptotically a complete 50(5) symmetry in the neighborhood of some critical point. 

In this paper, we show why symmetry-breaking effects due to the projection are asymptotically irrelevant in the 
neighborhood of a finite-temperature critical point. However, two kinds of symmetry-breaking effects tend to prevent 
50(5) symmetry from being restored asymptotically. One is related with the different mobilities of hole pairs and 
magnons (77 ^ 1 below), and the second one is due to the renormalization effects from quantum fluctuation at an 
intermediate length scale. The common tendency of these effects is to draw the system into a region of instability, whor e 
the two AF/normal(N) and SC/N transitions become first order before merging at the AF/SC/N triple pointc2lli3. 
However, when the first effect is large, quantum fluctuations tend to take the system back to the 50(5) point. 

This paper is organized as follows: In Sec. [n| we start from the projected 50(5) -symmetric model (allowing for a 
symmetry-breaking term 77 in the mobilities), and treat it by a slave-boson functional- integral approach in order to 
deal with the hard-core constraint. The important result is that at the AF-SC transition, the classical part of the 
action of the projected model preserves its 50(5) structure at n = I, despite the symmetry-breaking terms arising 
from the projection. These terms only appear in the quantum-mechanical (i. e., time derivative) part of the action. 
This fact gives a rigorous, justification for the much used semiclassical description of the high-T c materials via a 
50(5) -symmetric modeH'E3Eil despite of the presence of the large Hubbard gap. In Sec. [II, we derive the associate 
effective Ginzburg-Landau model by integrating out the momenta conjugate to the AF superspin variables£3. We 
study the properties of such model in the neighborhood of the AF/SC/N triple point and discuss the possibility of 
50(5) symmetry restoring at long wavelengths. 

In Sec. IV, we evaluate the corrections to the effective classical action due to the so far neglected quantum fluctua- 
tions. For small temperatures, these mainly affect the magnon-magnon scattering, thus breaking the 50(5) symmetry. 

Finally, in Sec. we draw our conclusions. Some details of the calculations are given in the appendices. 



II. THE MODEL 



We start from the effective bosonic model introduced in Refs. |2^JfJ, which describes low-energy bosonic excitations 
of "blocks" , (also referred to as "sites") labeled by the coordinate x, consisting of a rung in a I-D ladder or of a 2 x 2 
plaquette in a 2-D system. 

H = A s tt (x)t a (x) + A c A (x)U (x) 

X X 

- J s 2_j n a (x)n a (x') - J c 2J ni(x)m(x') (1) 

<xx'> <xx'> 

In this paper, we shall use similar conventions as in Ref. [l4|, where the indices a, b,.. are the 50(5) superspin 
indices and take the values 1,2,3,4,5 (in some cases, they might also include the "hole" index h), a, (3,.. — 2,3,4 
(corresponding to x,y,z) denote the spin indices and i,j = 1,5 denote the charge indices, and repeated indices are 
implicitly summed over. A boldface sign indicates the vector as a whole. Here, A s is the energy required to produce 
a magnon excitation, i. e. to replace a singlet with a triplet in a block, while A c is the energy required in order to 
produce a particle or hole pair. It is clear that A c is of the order of the Mott-Hubbard gap and, thus, A c ^> A s . J s 
and J c describe the hybridization of these excitations between nearest-neighbor sites and are related to their mobility. 
The Hamiltonian Eq. ([[]) acts on a "vacuum" |S1), which is a kind of "RVB" state consisting of a product state of 
half-filled singlet states |^(a;)) in each block. On the other hand, the five-fold states tj,(x)|£l(iz;)l-describe the triplet 
magnon states (for a — 2,3, 4), and the <i-wave hole and particle pair states on a block (a — 1, 5)E3. More specifically, 
one can define the charge eigenoperators t/j and t p as 

t x = —j=(t h + t p ) t 5 = —j=(t h - t p ) , (2) 

where t* h is the creation operator for a hole pair and t p is the creation operator for a particle pair. In Eq. (|l]), the 
n a play the role of the "displacement" coordinates of the local harmonic oscillators, while we will denote with p a the 
conjugate momenta, and we have the transformation to canonical variables: 
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t a = (n a + ipa) ■ (3) 

Due to their microscopic origin these bosonic states are hard-core bosons, in the sense that at most one boson can 
reside on each site. Mathematically, this is expressed by the condition 

ti(x)t a (x) < 1 . (4) 

The particle and hole density is controlled by the chemical potential /i, which couples to the Hamiltonian via a term 

H„ = -2^[tl(x)t p (x)-tl(x)t h (x)] . (5) 

X 

In the presence of this chemical potential term, the gap energy of the hole and particle pairs become A c + 2/x and 
A c — 2fi respectively. A (negative) chemical potential of the order of the charge gap A c /2 is needed to induce an 
AF-SC transition in this system. Near such a transition point, the gap energy of the hole pair A c + 2/z can be 
comparable to the (local) spin gap A s , while the gap towards a particle pair excitation is pushed up and becomes of 
the order of twice the charge gap. Since this is a very large energy scale, we can safely project this excitation out of 
the spectrum in the low-energy limit, by requiring that the condition 

*pO«OI*>=0 (6) 
is fulfilled at every site x. The new Hamiltonian takes the form 

H = A S ^2 tl(x)t a (x) + (A c + 2a*) ^ tU^hi*) ( 7 ) 

x,a X 

-J s ^2 n a (x)n a (x f ) 

<x,x'> ,a 



-j c /2 J2 (AixW) 



In Ref. II it was shown that the constraint Eq. (|^) can be enforced by introducing canonical commutation rules 
between the two variables m and ns, i. e. 

[m,n 5 ]=i/2, (8) 

and therefore we can identify \f%i\ with the "hole displacement" n>h and \[2n^ with its conjugate momentum ph- 
The 15*0(5) structure of the Hamiltonian becomes now clear if one introduces the superspin vector 

m a ee (r]n h ,n2,n 3 ,n4,rjp h ) , (9) 

where, for convenience, we have absorbed the different mobility for hole pairs and magnons rj = yj -^j- into the 

definition of the superspin. Carrying out the transformation to canonical variables Eq. ([}]), the Hamiltonian Eq. (0) 
now takes the simple form 

XX X 

~ J X! m a(x)m a (x') , (10) 



A c ee Ac "t 2M A s ee A s , and J = J s . (11) 



where we have further redefined 



The anisotropy in superspin space due to rj reflects now into the constraint, as we will see in Eq. ( |16| ) below. If 
one forgets for a moment the connection between coordinates m a and their conjugate momenta, the Hamiltonian 
Eq. ( |l0| ) becomes exactly 50(5) invariant under rotation of the superspin Eq. (^) at the AF-SC transition point 
A s = A c , which is reached by changing the chemical potential /i, i. e. at the AF-SC transition. If one further has 



3 



rj = 1, i. e. 2J S — J c , the constraint is invariant as well, and one apparently has a complete 50(5) symmetric model 
(cf. Ref. U). More specifically, one would like to 50(5) "rotate" just the m a coordinates, leaving the conjugate 
coordinates to the magnon part p a unrotated. This is possible, for example, in a classical ensemble, where, due to 
Liouville's theorem, expectation values are evaluated as J JJ. dpi dqi exp —H[pi, qi], and rotations of the qi only leaves 
the measure invariant. Of course, this does not hold for dynamics, which is affected by the relation between the two 
"superconducting" canonically conjugate components mi/77 and m^/rj, and between the AF components m a and their 
conjugate momenta p a . Thus, 50(5) symmetry is broken in dynamics, as pointed out in Ref. [l4|. Unfortunately, the 
relation between conjugate variables is also important in quantum-mechanical static averages, so that ground-state 
or finite-temperature averages are generally expected to break the symmetry when the full quantum problem is taken 
into account. 

In order to understand the nature of the symmetry-breaking terms, it is convenient to go over to a functional-integral 
representation of the partition function for the Hamiltonian Eq. (Q). The hard-core constraints can be conveniently 
taken care of by means of a slave-boson representation^, where the boson operator e(x) labeling "empty" sites is 
introduced. The detailed procedure is shown in Appendix |A|. After this transformation, the action takes the form 



S = Sqm + Scl , (12) 



where 



f 

Sqm = dr V 
Jo 



1 

-ip a (x,T) m a {x,T) ? m 5 (x,T) ihi(x,T) 



(13) 



(m Q indicates the time derivative of m a ) has the well-known form pq of the Feynman path integral, the i coming from 
the imaginary-time representation. Moreover, 

S C l = J dr ^p Q (a;, r) 2 + ^^m Q (a;,r) 2 + y^m.fi.r) 2 

X X X 

- ■/ ^ e(x, t) m a (x,T) e(x',T) m a {x' ,t)\ , (14) 

<xx'> 

where we have to replace 



/ p a (x,T) 2 m a (x,T) 2 m t (x,T) 2 
e{x,T) = il — — , (15) 
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which implicitly includes the condition 

p a (x,T) 2 m a (x,T) 2 mi(x, 



2r] 2 



< 1 , (16) 



and where W e have already carried out the transformation to canonical coordinates, Eq. (y), for the corresponding 
fieldsES. Eq. Jl4| ) is the correct classical limit of a projected, i. e. of the physical SO(5) model. Notice that the 
effects of the hard-core constraint is to introduce a renormalization of the boson hopping, and to bound the superspin 
magnitude, without, however, fixing its lengthen. Thus, the requirement that the superspin magnitude be unity should 
not be taken as a rigorous constraint of the SO(5) theory, at least not of the projected one (which is the physical one). 
On the other hand, one expects that in thcrhpmogcncous ordered phase this constraint might be a good assumption. 
A similar result has been shown by WegnerEZI, namely, that the orthogonality constraint in the exact SO(5) model is 
not a rigorous constraint, but it is favored at high temperature, as it maximizes the entropy. 

Eq. (|l^) clearly identifies the 5*0(5) -symmetry breaking terms. The classical action Scl is exactly SO(5) invariant 
at the AF-SC transition (A s = A c ) and for 77 = 1, while apparently incurable symmetry-breaking terms come from 
the time-derivative terms in Sqm- More specifically, with these values of the parameters, if one carries out an 
50(5) rotation within the superspin vector, Eq. (|9|), Scl remains invariant, while Sqm is changed. If quantum 
fluctuations are neglected, one can choose time-independent fields and set Eq. ([13]) to zero. In this case, any 
equilibrium expectation value is exactly 50(5) invariant. More specifically, let us take-a generic 50(5) rotation 
matrix R(n) = expin a r a parametrized by the vector n (r a are the 50(5) generatorsQ) , and f[m(x),p(x)} is a 
function of the superspin vector m(x), and, possibly, of p(a;). Then, the classical expectation values <>cl have the 
property 
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< f[*n(x),p(x)] >cl=< /[R ■ m(x),p(x)] > C l , 



(17) 



which is the requirement of SO(5) invariance. Notice that the p a should not be rotated, while in an exact 5*0(5) 
model they should. 

The question is: when is it justified to neglect the time dependence of the fields ?. This is allowed at moderately 
high temperatures, more precisely, at temperatures much larger than v/£ (in units of ks = h — 1), where £ is the 
correlation length and v is a typical velocity, in our case equal to J a, a being the lattice spacing. This means that 
neglecting Sqm is exactly justified when £ becomes infinite, i. e. in the neighborhood of a finite-temperature critical 
point, as a possible (finite temperature) multicritical point at which the AF/N and the SC/N transition lines merge 
into a first-order linea. Moreover, this critical point is indeed a good candidate for a possible asymptotic restoring 
of the complete SO(5) symmetry even in the presence, microscopically, of a projected SO(5) symmetry. This is 
very important as it would mean that the large-energy symmetry-breaking effect of the Mott-insulator gap would be 
exactly compensate d at . this critical point. This is analogous to the well-known situation for the antiferromagnetic 
spin-flop transitionar 8 lr ( Tr 9 rl 30 l where a system with uniaxial anisotropy restores SO(3) symmetry at the bicritical 
point. However, there are some important differences with respect to the spin-flop transition, as we will show in 
the next Sections. Moreover, notice that due to the symmetry breaking term, Eq. (fl3|), it is unlikely that SO(5) 
symmetry can be restored if the AF-SC transition is controlled by a quantum-critical point. Since we are interested 
in finite-temperature critical points, we will restrict to the case of three spatial dimensions D. 



III. EFFECTIVE GINZBURG-LANDAU ACTION 



In this Section, we study the action Eq. (|10|) in more detail. We first integrate out the momenta p a and obtain an 
effective action restricted to the superspin variables. For temperatures smaller than the singlet-triplet splitting A s , 
one can restrict to a Gaussian integration of the momenta, i. e., consider only quadratic terms in p a . Carrying out 
such an expansion, one obtains 

S C L = S pm + S m + 0(pi) , (18) 

where, leaving the r dependence implicit 

S m = [ d t [^^2m a (x) 2 + ^^2mi(x) 2 
J o 1 x z x 

- J J~] r(x) m a (x) r(x') m a (x')] , (19) 

<XX f > 



S pm = J dr^2 ^ A(x) ppix) 2 



(20) 



where we have defined 



9 T ( \ _ n " 

#)slt ^?g±-j-J^m.(x + d)r(x + d), (21) 



and the sum Y^d™ extends over nearest-neighbor sites. 

It is now convenient to reabsorb the m-dependent coefficient A(x) of the p 2 term into the definition of the momenta 
p. This is done in order to avoid the appearance of terms depending on the amplitude of the imaginary-time slice 
in the effective action. Furthermore, in order to avoid a m-dependent Jacobian due to the transformation, it is 
convenient to transform the m-coordinates in such a way that the Jacobian remains unity. The general procedure is 
illustrated in Appendix [b]. Up to second order in m 2 , the new m' coordinates are related with the old ones via 

3 7 

m a (x)=m' a (x)(l + — |m'(x)| 2 ). (23) 
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After this transformation, the integration of the p' a (x) 
QM action in the form 



p a (x) y/A(x) only affects Sqm, and one obtains a new 



S'r 



QM 



dr / dx 



2A S A(x) t] 2 



m${x) rhi(x) 



(24) 



where the transformation Eq. ( |23| ) should be inserted, and we have absorbed the unit cell volume V = a 3 in the 
definition of the fields by renaming m^/V — > m 2 



Thus, the total effective SO (5) action restricted to the superspin variables is given by Eq. (|19|) plus Eq. fl24|). The 
transformation Eq. (^) must still be carried out on the m variables, but, due to the fact that the coefficient tM- 
is small at the transition, this does not change the result significantly. On the other hand, it is important to take 
into account the effects of the hard-core constraint, which introduces the transformation Eq. (22), and, implicitly, 
the restriction of the superspin within a 5-dimensional hypersphere (or a ellipsoid, if 77 ^4-l)E3. Thus, S m ( Eq. (pj|)) 
gives the first effective classical functional microscopically derived from an 50(5) modeled, where the physics of the 
Mott insulating gap has been properly taken into account via the projection. This is the appropriate functional which 
should be used for physical predictions of the 50(5) theory, consistent with the gap. 

Close to the phase transitions, it is more convenient to derive a Ginzburg-Landau form for the action, obtained, as 
usual, by expanding in powers of the field m and keeping only lowest-order gradient terms. After inserting Eq. (|22" 
and dropping the prime indices in the fields m, we obtain 



dr / dx 



\j^rn a (x) 2 + ^■m i (x) 2 + ^(Vm a (x)) 2 
%fcma(tf\ +^fl>( a 9 ) +x(5>«(*) 2 ) (E^(*) 2 )}> 



(25) 



with the parameters 



T s 1 c. , 1 c 

2 ~ ^ 2 


DJ) 


p J a 2 




2 2 




8 2 V 


3 r s , 
7 A/ 


u c V J D 
~8 ~ ~2~V 4 


3 Tc . 


7 A s 


U c + U s 
Ucs - 2 




(p5|) have been 


analyzed 



(26) 



The critical properties of the model Eq. (]25|) have been analyzed in several works§S00S@@. Its phase 
diagram is determined by two relevant parameters, the first one r s — r c oc A s — A c controls the transition between 
the AF and the SC phases, while the other ~ min(r s , r c ) controls the second-order transition between the appropriate 
ordered (AF or SC) andihe normal phase. At the transition point r s ~ r c ~ 0, there are two competing fixed point 
controlling the transitionEJ, the Heisenberg bicritical fixed point (in this specific case, the SO(5) fixed point), and the 
biconical tetracritical fixed point. According to the e-expansion, the latter fixed point turns out to be the stable one 
for n > n c ~ 4 — 0(e). This means that, in general, the model Eq. (|2f]), which has n = 5, is expected to flow to this 
latter fixed point and not to the 50(5) -symmetric one for u s ^ u c ^ u cs - On the other hand, since n = 5 is not very 
far away from n c , the stable biconical fixed point is approximately 50(5) invariant with symmetry- breaking terms of 
the order of 20%. Moreover, therais-a plane in the u s , u c , u cs space, given by the condition u 2 s = u c u Sl from which the 
system flows to the 50(5) pointEHrEJ. This is due to the fact that a scale transformation of, say, the SC components 
m 2 — » m 2 u s /u c of the order parameter would yield again an 50(5) -symmetric interaction of the form u|m| 4 . The 
asymmetry would then be transfered into different susceptibilities p s ,p c for the AF and for the SC order parameters. 
However, it has been shown in Refs. [33|,^9] that the different in the susceptibilities is an irrelevant parameter. 
In our case, we have 



Au 2 



> 



(27) 



G 



which means that the 50(5) symmetric fixed point is never reached, except when the equal sign holds, i. e., when 
rj 7^ 1 (at the transition r s — r c ). On the other hand, we expect on physical grounds the mobility of the hole pairs 
to be smaller than that of the magnons, and, thus, 77 to be smaller than 1. Unfortunately, for the case Eq. (|27|), the 
couplings flow away into a region of instability. The common interpretation is that the AF/N and SC/N transitions 
becornfij-first order as well (fluctuation-induced first-order transition), at least close enough to the AF/SC/N triple 
pointEM. 

This fact seems in contrast with the apparent observation of bicritical behavior with 50(5) .critical exponents in 
the organic superconductor k-(BEDT-TTF) 2 X (see Refs. |34] , j35|) , by Murakami and NagaosaHEEl. There may be 
several ways to understand this. rQne possibility is that other effects not considered here, such as, e. g., Coulomb 
interactions, fermionic excitationsEj, or quantum effects, as discussed in Sec. IV, counterbalance this effect and draw 
the system back to the domain of attraction of the biconical fixed point. As discussed above, the differences between 
the 50(5) and the biconical fixed point are only about 20%, so that they might be not observable experimentally. 
Alternatively, since the flow would <j*ft*j the SO (5) plane, it could produce SO (5) exponents at intermediate length 
scales. On the other hand, Hu et al.E§e3, observe a coexistence region of AF and SC for the 50(5) -anisotropic case, 
which could be possibly identified with the biconical phase. Their result could be due to the fact that they consider 
a different c-axis anisotropy (x in footnote Eoh , for the AF and for the SC variables. 



IV. QUANTUM CORRECTIONS 

Even when considering a classical (i. e., finite-temperature) critical point, the quantum-mechanical symmetry- 
breaking terms Sqm although irrelevant in the RG sense, contribute to the RG flow up to a certain length scale of 
the order of v/T. Since Sqm breaks the 50(5) symmetry, it is expected, during this initial renormalization process, 
to introduce symmetry-breaking terms in Scl- Therefore, even when Scl is 50(5) symmetric at the microscopic 
scale, the renormalized Scl at the scale £ ~ v/T will probably break the symmetry. In this Section, we evaluate 
these symmetry-breaking terms originating from Sqm, or, more precisely, from the time dependence of the fields. 

In order to evaluate these effects, we separate the fields into their static and dynamic parts, and integrate out 
the latter. Since we are working at finite temperature, we have to integrate out the components of the fields with 
Matsubara frequencies uj n = 2imT with n 0. In order to obtain an analytic expression for these corrections, we 
restrict to one-loop contributions and take just the leading low-temperature terms. 

We first diagonalize the non-interacting (quadratic) part of the action Eq. ((2J) plus Eq. ( ^oj ) by Fourier transform. 
We can neglect the corrections to S'q M due to the transformation to the primed variables Eq. (p3[), as it introduces 
irrelevant quartic time-derivative terms. In Fourier space, the action takes the usual form 



Sqm + S'cl = H rn a (-k) [G(fc)- 1 ] ab m b {k) + ±£ » 

z ^k °^k 1 ,k 2 M 



m a (ki) m a (k 2 ) u ab m b (k 3 ) m b (-ki - k 2 , -fc 3 ) , (28) 

where we have introduced the shorthand notation k = (k, u>), and ^£ fc = jj J2 U J A > with A ~ I/a a short-distance 
cutoff for k. In Eq. (p8|), the nonzero elements of the (non interacting) Green's functions read 

G ^ = JSA h2 > (29) 



and 



g(t)l , 1 = G ( fc ) M > ^ +4 , (30) 



G(k) 5A = -G(fc) 1 , 5 = -r^-r— t , (31) 



(pk 2 + r c ) 2 + 



and the interaction parameters are u a _p = u s , u.^j = u c , and Ui :Ct = u cs . 

At one loop, integration of the w / fields only changes the parameters r, and u, similarly to conventional field 
theory. In the T — > limit, the change of the former is finite, and merely shifts the transition point. On the other 
hand, the changes Su a b in the interaction parameters u a b grow logarithmically with decreasing temperature at the 
critical point. We will, thus, restrict to evaluation of these corrections. These are given by the sum of the usual "loop" 
diagrams, which give 
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Su ab — - - ^2 u ac u cb I cc - 2 u 2 ab I ab - u ab (J aa u aa + I bb u bb ) , (32) 



where the integrals I ab are given by 



G(k) aa G(-k) bb . (33) 



In Eq. ( |32] ) and Eq. (|33"|), we have neglected contributions from nondiagonal parts of Green's functions Eq. d3l|), 
as they only give finite contributions to integrals of the form Eq. ( |33| ) in the low- temperature limit. The same holds 
for integrals containing at least one Green's function of the. superconducting fields Eq. (^0|). This is due to the fact 
that for these fields the (bare) dynamical critical exponentcd'Cj z is equal to 2, and it does not produce divergences in 
D = 3. This in turns occurs because the two components of the SC order parameter are canonically conjugate, while 
the AF ones have independent massive ones. Therefore, we will consider only the divergent contribution 




A s AV^7 

~ 2ttT ' (34) 



where we have assumed that we lie outside of the region of influence of the quantum critical fixed point, i. e. 



T ^> y/r c / s A s . Replacing Eq. ( |34| ) in Eq. (|32j), we obtain for the leading contributions 

Su s = ~ u 2 s I„ (35) 
Su c = -| u 2 cs I s 

c _ _5 

As expected, quantum fluctuations draw the system away from the 50(5) -invariant point even in the case where 
rj = 1. This can be seen by adding these corrections to an initially 50(5) -invariant system with u c — u s — u cs = u. 
At the lowest order in the u a , the renormalized parameters u' a = u a + 5u a obey the relation 

Am' 2 = u' c 2 s - u' s u' c = 2 I s u 3 > , (36) 

i. e., as in the case of r\ ^ 1, Eq. (p7|), the system is drawn into the instability region where a fluctuation- induced 
first-order transition is expected. This indicates that quantum fluctuations and anisotropy 77 ^ 1 cooperate in the 
same direction and draw the system into the instability region, where no finite fixed point is expected. However, for 
the case where the u a are different, one obtains 

Au 2 ~ Au 2 = (7u 2 cs - 11 u c u s ) I s . (37) 



Further inserting the values of the u a from Eq. (|26j) with rj ^ 1 (we fix ourselves at the triple point r s = r c ), Eq. ( |37 
becomes negative for 77 < x c , or 77 > l/x c with x c ~ 0.498. Therefore, for large difference in the mobilities 77, quantum 
fluctuations tend to shift the renormalized parameters back towards the domain of attraction of the biconical and of 
the 50(5) fixed point. 



V. CONCLUSIONS 



In conclusion, we have analyzed the properties of a projected 50(5) model which takes into account the high-energy 
physics of the Mott-insulating gap. As already pointed out in Rcf. |l4|, the chemical potential can always be shifted to 
the AF-SC transition point in order to cancel the symmetry-breaking terms produced by the gap in the classical part 
of the action. On the other hand, symmetry-breaking terms due to the projection show up in the quantum-mechanical 
part of the action, as a conjugacy relation between the superconducting components of the superspin vector. A further 
source of symmetry breaking is due to the different mobility of the hole pairs and of the magnons parametrized by 

Close to the AF/SC/N finite-temperature multicritical point, the quantum effects due to the projection are irrel- 
evant, although subleading symmetry-breaking corrections appear at intermediate length scales. When considered 
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separately, these symmetry-breaking effects both draw the RG flow into a region of instability with first order transi- 
tions and no SO(5) symmetry. On the other hand, for strong anisotropies r\ < 0.5, quantum corrections partly cancel 
the symmetry-breaking effects. 

There are possibly other effects, such as Coulomb interaction, or fermionic excitations, which can possibly take 
the system back into the domain of attraction of the biconical fixed point, where SO(5) symmetry is only broken 
by ~ 20%. Notice that, since the order parameter must be rescaled in order to reach this fixed point, the (possibly 
approximate) 50(5) symmetry reached at this critical point is renormalized, in the sense of Rcf. [l^. This means, for 
example, that the S'0(5X-picture would be consistent with different absolute magnitudes of the SC and AF gaps, as 
observed experimentallytJ. 
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APPENDIX A: EXACT SLAVE-BOSON TREATMENT OF THE CONSTRAINT 



The hard-core constraint Eq. ([|) becomes (after projecting out the electron pairs) 

Q{x) = tl{x)t a {x) + t{{x)t h {x) + e\x)e(x) -1 = 0. 
The "physical" bosonic operators are then obtained as usual by the replacement 

t a (x) -> t a {x)e\x) , 



(Al) 



(A2) 



(including a = h) s o that the constraint is now conserved by the Hamiltonian. Within the functional integral, the 
constraint Eq. (Al) can be enforced as usual by adding a "Lagrange multiplier" term i^2 x X(x)Q{x) and integrating 
over all \(x). The partition function can thus be written in terms of an integral over bosonic fields 



Z = J Vil Vt a Vt{ Vt h De f Ve dXexp-S' 



with the action 



S' = jf dr {J2 r) (J^ + iA(a:)) t a (x, r) + t{{x, r) ( 



d_ 



iX(x) th(x,r) 



(A3) 



(A4) 



+e^(x, r) ( — — h iX(x) ] e(x, r) — i\(x) 



0t 



H(t 



>} 



where H(t) is obtained by replacing Eq. ( A2) in Eq. (0) and by replacing all bosonic operators with the corresponding 
fields at the imaginary time r (since the Hamiltonian is already normal ordered). In principle, one-should take a 
discretization of the time variable and consider the continuum lim it o nly at the end of the calculatiorJi3c30 Notice 
that the integration of A would not give a constraint like Eq. (Al) for the bosonic fields at all imaginary times. 
Nevertheless, one can proceed in the usual way by carrying out the gauge transformation 



e(x,T) = e(x,T)e ze{x ' T) 
t a {x,T)=t a {x,T)e ie( - x ^ 
\{x) = X(x, t) — 9(x, t) , 



(A5) 



where e(x, r) = \e(x, r)|. In this way, we can restrict to real values of the boson filed e and absorb the time dependence 
of its phase into a (now) time-dependent A. Integration over X(x, r) now leads to the enforcement of the constraint 
via the 5 function (for simplicity, we drop the bar everywhere) 



JJS [\t a (x,T)\ 2 + \t h (x,T)\ 2 +e(x,r) 2 - l] 



(A6) 



at all imaginary times. Integration over e(x, r) allows one to replace it everywhere in the Hamiltonian, leading to the 
new action Eq. @ with Eq. (0). 
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APPENDIX B: INTEGRATION OF THE MOMENTA 

The p-dependent part of the action has the general form (Cf. Eq. ( [j~3"| ) and Eq. (pOj)) 

S p = f dr [ dx A A(\m(x)\ 2 ) Pa (x) 2 - ip a (x) B{x) , (Bl) 



where A is a function of the superspin's magnitude squared (for simplicity, we neglect gradient terms). In order to 
absorb the coefficient A, we define new momentum variables 



p> a (x)=p a (xWA(\m(x)\2). (B2) 

However, since we don't want to produce an m-dependent Jacobian, we carry out a similar transformation for the m 
variables as 

m' a (x) = m a (x)g[\m(x)\ 2 } , (B3) 
where g is chosen in order to have a Jacobian equal to 1. This requirement gives the differential equation 

A(|m| 2 ) 3 / 2 [ 5 (|m| 2 )" - 2|m| 2 3 (|m| 2 )- V(|m| 2 )] = 1 , (B4) 
n (= 5) being the number of components of the superspin m. The solution of this equation is 

(Vfg(r)T = ~j r^^Air)^' 2 dr , (B5) 



where r = |m| . Upon restricting to the lowest order of Eq. (21), A(r) = 1 + TK~ r + )> we obtain 



3 / 

g[\m(x)\ 2 ]^(l-— \m(x)\ 2 ), (B6) 



and its inverse Eq. (E3| 
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